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Orbital magnetism is studied for graphene flakes with various shapes and edge configurations using 
the tight-binding approximation. In the low-temperature regime where the thermal energy is much 
smaller than to the energy level spacing, the susceptibility rapidly changes between diamagnetism 
and paramagnetism as a function of Fermi energy, in accordance with the energy level structure. 
The susceptibility at charge neutral point is generally larger in armchair flake than in zigzag flake, 
and larger in hexagonal flake than in triangular flake. As the temperature increases, the discrete 
structures due to the quantum confinement are all gone, and the susceptibility approximates the 
bulk limit independently of the atomic configuration. The diamagnetic current circulates entirely on 
the graphene flake at zero temperature, while in increasing temperature it is localized near the edge 
with the characteristic depth proportional to 1/T. We predict that the diamagnetism of graphene 
can be observed using the alignment of graphene flakes in a feasible range of magnetic field. 



I. INTRODUCTION 

The recent developments in fabrication of graphene- 
based systems realized a variety of graphene nanos- 
tructures, such as graphene ribbons^ and graphene 
flakes j£r— The electronic band structure in these sys- 
tems crucially depends on the shape and also on the 
edge termination,— lii - — giving physical properties dis- 
tinct from those in bulk graphene. So far, a number of 
theoretical researches have been devoted to understand- 
ing the electronic properties of graphene ribbons 1 *^— and 
graphene flakes 1 ^ - — with various atomic configurations. 



(a) Hexagonal Armchair (b) Trigonal Armchair 




(c) Hexagonal Zigzag (d) Trigonal Zigzag 




FIG. 1: Atomic structures of (a) hexagonal armchair, (b) trig- 
onal armchair, (c) hexagonal zigzag, and (d) trigonal zigzag 
graphene flakes. 



The purpose of this paper is to study the orbital mag- 
netism of graphene flakes. Experimentally, the mag- 
netic property of graphene-based materials was investi- 
gated for bulk graphite— - — , nanographite^fi, and exfoli- 
ated graphene nano crystals^. There the susceptibility 
always contains a strong diamagnetic background due 
to the orbital effect, whereas it is also contributed by 
the spin paramagnetism, 31 and in some cases the sponta- 
neous spin magnetic ordering 2 ^ 30 ' 32 which can be caused 
by the zero-energy edge state o 13 ' 33 ! 34 and atomic defects. 
In any case, correct understanding of the orbital suscepti- 
bility of finite graphene systems is important to describe 
the overall magnetic property in realistic graphene sys- 
tems. 

Graphene has unusual electronic band structure char- 
acterized by the massless Dirac spectrum^r— and ac- 
cordingly, the orbital magnetism is significantly different 
from the conventional Landau diamagnetisnij 38 i 44 ~ — The 
orbital susceptibility of bulk graphene diverges when the 
Fermi energy resides at Dirac point, but vanishes inside 
the conduction or the valence band. Finite-size effect on 
this singular diamagnetism has been theoretically studied 
for carbon nanotubes^— and graphene ribbons i 13 i 59 In 
our previous work,— particularly, we found that the sus- 
ceptibility of graphene ribbon behaves in a complicated 
manner as a function of Fermi energy, reflecting the sub- 
band quantization imposed by the spacial confinement. 

In this paper, we consider the orbital diamagnetism 
of lower dimensional systems — graphene flakes as illus- 
trated in Fig. [1] For each case we calculate the orbital 
magnetic susceptibility and the diamagnetic electric cur- 
rent distribution using the tight-binding model. We find 
characteristic properties peculiar to each different case, 
and also general tendencies independent of the configura- 
tion. We also predict that the diamagnetism of graphene 
can be observed using the alignment of graphene flakes 
dissolved in a solvent under a magnetic field. The paper 
is organized as follows. In Sec. HH we introduce tight- 
binding Hamiltonian and the formulas to describe the 
orbital magnetic effect. We calculate the magnetic sus- 
ceptibility and the diamagnetic electric current distribu- 
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tion for graphene flakes in Sec. Mil and IIV1 respectively. 
We make a quantitative comparison between the orbital 
magnetism and the spin magnetism in Sec. |Vl We argue 
the magnetic-field alignment effect in Sec. lVII and present 
a brief conclusion in Sec. IVIII 



For the later references, let us review the orbital mag- 
netism of the bulk graphene. The low-energy physics of 
graphene electrons can be effectively described by the 
continuum massless Dirac Hamiltonian 3 -^— and the or- 
bital susceptibility is calculated for this model a a 38 i 48 ' 53 



II. FORMULATIONS 

Graphene is composed of a honeycomb lattice of car- 
bon atoms, where a unit cell contains A and B sublat- 
tices. We consider four different atomic configurations 
of graphene flakes as shown in Fig. [TJ which are charac- 
terized by hexagonal or trigonal shape and by armchair 
or zigzag edge termination. For each case, we range the 
system size from a few nm to a few tens of nm. We de- 
scribe the motion of graphene electrons using the nearest- 
neighbor tight-binding model for p z atomic orbitals. The 
Hamiltonian is written as 



H = 



-7o 



4, 



(1) 



(n,m) 



where —70 is the transfer integral, c\ is the creation oper- 
ator at the site n, and (n, m) represents summation over 
all nearest-neighbor sites. The parameter 70 was exper- 
imentally estimated in the bulk graphite as 70 ~ 3eV. 60 
The system is under a uniform magnetic held B perpen- 
dicular to the graphene plane, which is incorporated by 
the Peierls phase 4> n m, 



t>nm = —r I d£-A, 
cn 



(2) 



where A(r) is the vector potential giving the magnetic 
field by B = V x A. 

For each single graphene flake, we diagonalize Hamil- 
tonian Eq. (JTJ and obtain a set of eigenenergies £j. The 
thermodynamical potential at temperature T and chem- 
ical potential \i is written as 



(3) 



The magnetic susceptibility per unit area is given by 



X 



S\dB*)^ T 



(4) 



B=0 



where S is the area of the system. To calculate this, 
we obtain the eigenenergies at zero magnetic field and a 
small finite magnetic field, and numerically calculate the 
derivative of the thermodynamic potential. 

The electric current from the site m to n is calculated 
by an operator, 



J n 



* c ri c m 



h.C. 



(5) 



We obtain the expectation value of J mn for each bond 
using the eigenstates at a sufficiently weak magnetic field, 
where the current amplitude behaves linearly to B. 
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24ttc 2 k B T cosh 2 \p,/(2k B T)] 



(6) 



where g v — g s = 2 are the valley and spin degeneracies, 
respectively, v is the band velocity related to the transfer 
integral by hv = \/3a7o/2, and a w 0.246nm is the lattice 
constant of graphene. At T = 0, Eq. © becomes a delta 
function, 



Xes(^,T = 0) = -g v g. 



e v 
6nc 2 



(7) 



III. MAGNETIC SUSCEPTIBILITY 



Fig. [5] shows the susceptibility against the chemical 
potential for four types of the graphene flakes with several 
different temperatures. The areas of the flakes are taken 
to be nearly equal to S « (23.5nm) 2 , which includes 
1.1 x 10 4 of hexagonal unit cells. The horizontal and 
vertical axes are scaled by 



£0 



Xo 



hv 

g v g s e 2 v 2 
6ttc 2 Sq 



(8) 
(9) 



respectively. £0 represents the energy scale in the Dirac 
cone associated with the length scale y/S. We also cal- 
culated the susceptibility for different system sizes and 
found that for each of four types, the susceptibility and 
the level structure plotted in this scale becomes almost 
universal as long as S> a. This is naturally expected 
from the fact that the low-energy physics are well de- 
scribed by the Dirac Hamiltonian. 

Upper figure in each panel presents the energy level 
structure at B = 0, where dashed (black) lines are non- 
degenerate levels, and solid (red) lines are two-fold degen- 
erate levels. In the low temperature regime, k B T <C £0, 
we observe that the susceptibility abruptly changes at 
every single energy level, and in particular, it exhibits 
sharp spikes toward the paramagnetic direction (down- 
ward in the figure) at two-fold degenerate levels. This 
is because the degenerate states, having opposite mag- 
netic moments, split linearly in magnetic field just like 
spin Zeeman splitting, and induce paramagnetism in an 
analogous way to spin Pauli paramagnetism. The con- 
tribution to the orbital susceptibility (per area) from the 
degenerate states at Eq is written as 



2to" 



6(11- E ), 



(10) 



where ±m is the magnetic moments of the doublet. The 
typical magnitude of m is shown to be y/Sev/c, which 
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FIG. 2: Orbital magnetic susceptibility as a function of /i in (a) hexagonal armchair, (b) trigonal armchair, (c) hexagonal 
zigzag, and (d) trigonal zigzag graphene flakes with the size of \fS ~ 23.5nm. In each figure, the upper panel presents the 
energy spectrum, where dashed (black) lines represent non-degenerate levels, and solid (red) lines two-fold degenerate levels. 



is the only magnetic-moment scale in the massless Dirac 
system. 

The major difference between armchair flakes and 
zigzag flakes comes from the existence of the zero-energy 
edge states peculiar to the zigzag edge.— ~— In the trian- 
gular zigzag flake, Fig.[2Jd), there are a number of energy 
levels exactly at zero energy^ of which wavefunctions are 
shown to be localized at the edge, and the degeneracy is 
the order of ~ \/~S/a. Remarkably, the susceptibility in 
the low temperature regime is completely flat at these 



levels, meaning that the edge states have absolutely no 
contribution to the orbital magnetism. This is simply be- 
cause the edge states are locked to zero energy even in the 
presence of magnetic field, and never participate in the 
total energy change. In the hexagonal zigzag flake, Fig. 
HJc), on the other hand, the edge levels slightly shift from 
zero energy, leading to some small contributions to the 
orbital susceptibility. The energy shift arises because the 
edge states on neighboring sides of the hexagon always 
reside at different sublattices, and they are hybridized 
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FIG. 3: Extended plot of the susceptibility curves in Fig. [2] 
over the whole band region, for the four types of graphene 
flakes at ksT/eo = 2.22. The energy axis is now scaled 
by absolute unit 70, in which the temperature amounts to 
feT/70 = 0.02. Inset shows the detail of the central peak. 
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by some finite matrix element including 70. Neverthe- 
less, the edgestates do not play a significant role in the 
overall behavior of the orbital magnetism. 

As the temperature increases, the spikes and steps 
in the susceptibility are smeared out into an oscilla- 
tory curve. The oscillation eventually disappears in 
hsT/so >^ 1, leaving a single diamagnetic peak at 
the Dirac point, which corresponds to the thermally- 
broadened delta- function in the bulk limit, Eq. ([5]). In 
Fig. [3l we present an extended plot of the susceptibil- 
ity curves at k^T/e^ = 2.22 over the whole band region, 
for different types of graphene flakes with y/S rs 23.5nm. 
The energy axis is now scaled by absolute unit 70. We 
see that the finite-size effect almost vanishes in this tem- 
perature regime, giving a universal curve independent of 
the atomic configuration. The curves still slightly differ 
near the central peak, because the level spacing around 
the Dirac point, which is about ~ Eq, is not completely 
negligible compared to k^T at this particular system size. 
This small variance would vanish in larger flakes which 
satisfy > £o- 

The curve is characterized by the strong diamagnetic 
peak at the Dirac point and some smaller structures off 
the Dirac point. The contribution from the lower-half 
spectrum adds up to a paramagnetic offset to the central 
diamagnetic peak. Namely, the susceptibility near /z = 
is approximately written as 
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FIG. 4: (a) x(A t = 0; T) of hexagonal armchair flakes with 
several different sizes, plotted in the absolute units 70 and 
g v g 3 e 2 v 2 /(6TTc 2 'y ). (b) x(p = 0;T) - Xpara of the same sys- 
tems, plotted in the relative units £0 and \o- (c) Plot simi- 
lar to (b) for different types of graphene flakes with the size 
23.5nm. 



where Xeff is given by Eq. ©, and 



Xpa 



0.37 x 



g v g s e 2 v 2 

67TC 2 7o 



(12) 



(11) 



The offset Xpara is much smaller than the height of the 
central peak ps g v g s e 2 v 2 / ' {2A-Kc 2 k-QT) 1 since k^T is usu- 
ally much smaller than 70. Outside the Dirac cone region, 
we see tiny Landau diamagnetism in the quadratic band 
bottom at \i = ±370, and paramagnetism around the van 
Hove singularity at \i = ±701^— 
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To analyze the size dependence quantitatively, we plot 
x(/i = 0; T) of hexagonal armchair flakes with different 
sizes in Fig.@|a) and (b). The panels (a) and (b) present 
the same information but in different fashions: (a) shows 
X in the absolute units 70 and g v g s e 2 v 2 / (67rc 2 7o) for hor- 
izontal and vertical axes, respectively, while (b) shows 
X ~ Xpara the contribution from the Dirac cone, with 
relative units £o and %o depending on the system size. 
In (b), we see that the curve converges to a single uni- 
versal curve as the size increases, indicating that the 
physics there is well described by effective Dirac equa- 
tion. The susceptibility approaches the bulk limit Xeff 
in the high temperature region k^T 3> £o, whereas in 
k^T <^ eq it deviates from Xoff and reaches some finite 
maximum value. When we consider the susceptibility 
of a single graphene flake, xS, at a fixed absolute tem- 
perature, it scales in proportion to xo^ oc S 3 ^ 2 in the 
low-temperature regime k^T <^ Sq, while it is just pro- 
portional to S in the high-temperature regime k^,T 3> £o 
where x is equal to the bulk limit. 

The detail of the universal curve in Fig. HJb) depends 
on the flake shape and the edge configuration. In Fig. 
IDJc), we present plots similar to Fig.UJb) for four different 
types of graphene flakes with the same size VS« 23. onm, 
which is sufficiently large to achieve the universal limit. 
In low temperatures, the susceptibility tends to be larger 
in an armchair flake than in a zigzag flake, and larger 
in a hexagonal flake than a trigonal flake. In the high 
temperature region, on the other hand, all the curves 
approaches the same bulk limit. A similar edge depen- 
dence was previously found in graphene ribbons, where 
armchair ribbons generally exhibit larger diamagnetism 
than zigzag ribbons . 13 ! 59 



IV. DIAMAGNETIC CURRENT 
DISTRIBUTION 

Fig. [5] shows the diamagnetic current distribution in- 
duced by the magnetic field in the four types of graphene 
flakes of the size VS ^ 23. 5nm at several different tem- 
peratures. To visualize the global current circulation, 
we illustrate continuous flux lines obtained by smoothing 
the original discrete current J mn on each bond, which is 
shown in the left inset. Specifically, we find a certain 
potential function ^ which satisfies J = e z x V^, and 
obtain the equi-potential lines of VP as the current flux 
lines. At zero temperature, the flux circulates entirely 
on the system reflecting the absence of the characteris- 
tic wave length in graphene. As temperature becomes 
higher, it is going to be localized near the edge. The cur- 
rent circulation of zigzag and armchair graphene flakes 
are globally similar, but the flux lines of armchair flakes 
exhibit some roughness while it is not observed in zigzag 
flakes. This actually corresponds to the atomic-scale cur- 
rent circulation in the Kekule pattern seen in the original 
current map (left inset), 13 which is caused by the inter- 
valley (between K and K 1 ) hybridization peculiar to the 



armchair edge. 

Fig. |5] shows the detailed plots of the electric current 
as a function of position from the boundary to the center, 
calculated for (a) the zigzag and (b) armchair flakes. The 
position is labeled by the bond index defined in the inset, 
and A and B (B') correspond to the edge and the center 
of triangle (hexagon) , respectively, which are depicted in 
Fig. [SI C )- The current distribution is more localized to 
the edge when T becomes higher, and the typical depth 
of the edge current is characterized by 

hv 

Acdgc = 2^T' (13) 

in accordance with the result for graphene ribbons^ 

The current distribution in the atomic scale strongly 
depends on the edge type. We can show that, how- 
ever, the integrated edge current approximates cXeffB 
independently of the edge type, in the high temperature 
regime fceT/eo > 2. This is consistent with the fact 
that the orbital susceptibility is then given by the bulk 
limit regardless of the atomic configuration. When com- 
paring hexagonal and triangular flakes of the same edge 
type, we see that the curves are almost completely equiv- 
alent in k^T/eQ > 2. This suggests that the edge current 
distribution in high temperature is solely determined by 
the local edge configuration, independently of the global 
shape. 

V. COMPARISON TO SPIN PARAMAGNETISM 

The orbital magnetism always competes with the spin 
paramagnetism which has been neglected so far. When 
we include spin Zeeman splitting, each spin-less energy 
level at Eq acquires the Pauli paramagnetism 

XPauii = ^(|^ B ) 2 2^ M -£; ), (14) 

where g ~ 2 is the g factor for a graphene electron, 
fiB = eh/(2nioc) is the Bohr magneton and too is the 
bare electron mass. This is similar to the orbital contri- 
bution of Eq. (TTUj) for doubly degenerate levels, while the 
orbital magnetic moment to there is now replaced with 
(?/ib/2. In the flakes of S > (lnm) 2 , /zb is much smaller 
than the typical magnitude of to, which is rs ySev/c, 
suggesting that the Pauli paramagnetic effect is typically 
much smaller than the orbital effect. This in contrast 
to conventional electron systems where orbital magnetic 
moment and spin magnetic moment are both of the order 
ofn B M 

In a zigzag graphene flake, the highly-degenerate edge 
states at zero energy give exceptionally large Pauli para- 
magnetism. The contribution is written as 

XPauii = ^(f/x B ) 2 25(A), (15) 

where iV c dgc(~ \/~S/a) is the number of edge states per 
spin. In the low-temperature regime such that k^T <C £q, 



6 



(a) Zigzag 

k B T/eo=0 k B T/eo=1.11 k B T/eo=2.22 




k B T/so=0 k B T/eo= 1 . 1 1 k B T/eo=2.22 




k B T/eo=0 k B T/eo= 1 . 1 1 k B T/eo=2.22 




FIG. 5: Diamagnetic current distribution in different types of graphene flakes of the same size %/S « 23.5nm at several different 
temperatures. Continuous flux lines are obtained by smoothing the original discrete current on each bond, which is shown in 
the left. 

this is dominant over the orbital effect near zero energy, it should be compared to the bulk orbital susceptibility 
since the orbital susceptibility does not diverge at edge Xcffi Eq. ([5]). The ratio between two opposite components 
states as already shown. In high-temperature regime 
kftT ^> £q, the delta- function is thermally broadened and 
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FIG. 7: Angle distribution of hexagonal armchair flakes in 
magnetic fields at T — 300K. 



FIG. 6: Electric current as a function of position from the 
boundary to the center, in (a) the zigzag and (b) armchair 
flake with y/S w 23.5nm.. The position is labeled by the 
bond index defined in the inset, and A and B (B') represent 
the border and the center of triangle (hexagon), respectively, 
which are specified in (c). Horizontal arrows indicate A c d gc 
for k B T/e = 1.11, 2.22 and 4.44. 



graphene-based materials^ -30 ' 32 The origin of the spon- 
taneous magnetism is still under debate, while it is sup- 
posed to be caused by the atomic defects, grain bound- 
aries, and highly-degenerate edge statesj 13 ' 33 i 34 
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so that the Pauli paramagnetism is negligible in a large 
flake with 3> fl- 
it should be noted that graphene flakes may have lat- 
tice vacancies and/or adatoms depending on the experi- 
mental condition, and the impurity levels given by these 
defects contribute to additional Pauli paramagnetism. 
Moreover, we remark that several experimental studies 
reported the evidence of ferromagnetic spin ordering in 



VI. MAGNETIC FIELD ALIGNMENT OF 
GRAPHENE FLAKES 



The diamagnetism of graphene can be possibly ob- 
served using the magnetic-field alignment of graphene 
nanoflakes dissolved in a solvent, similarly to the experi- 
ments for the carbon nanotube^ In a magnetic field, the 
graphene flakes tends to be oriented parallel to the field 
direction, because the field component penetrating the 
graphene plane raises the total energy due to the dia- 
magnetism. If we assume that the graphene flakes are 
planer and rigid, the condition to achieve the alignment 
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is roughly estimated as 

\xB 2 S > k B T. (17) 

For the graphene flakes VS « 23.5nm at T = 300K, for 
example, the required field becomes B > 9T. 

We calculate the angle distribution of graphene flakes 
with various sizes using the Maxwell-Boltzmann statis- 
tics. In the thermal equilibrium, the probability that the 
normal of the graphene plane is inclined from the mag- 
netic field by 9 to 9 + d9 is written as P(cos9)d(coa9), 
where 

expj-^(cos^ 

f\ exp[-/3U (cos 0)]d(cos 9) 

with U(cos9) = -{l/2) X SB 2 cos 2 9 and /3 = l/{k B T). 
Fig. [7] plots the distribution function P(cos#) calcu- 
lated for hexagonal armchair flakes with several sizes at 
T = 300K, using x in Fig. |4|a). We see that the align- 
ment occurs more strongly in larger flakes, because the 
magnetization of a single flake, x^B, is greater for larger 
S. Note that it is not only due to a linear factor S 1 , but 
also because x increases in larger S as shown in Fig.|4][a). 

VII. CONCLUSION 

We have studied the orbital diamagnetism of the 
graphene flakes with various shapes and edge configura- 
tions using the tight-binding approximation. We found 
that the behavior is significantly different depending on 
the relative magnitude of the thermal broadening energy 
k B T to the typical energy level spacing £o = hv/V~S. In 
the low-temperature regime where k B T <^ Sq, the suscep- 
tibility as a function of Fermi energy rapidly changes be- 
tween diamagnetism and paramagnetism in accordance 



with the level structure depending on the specific atomic 
structure of the flake. The susceptibility at the zero 
Fermi energy is found to be generally larger in armchair 
flakes than in zigzag flakes, and larger in hexagonal flakes 
than trigonal flakes. In the high-temperature regime 
k B T > 2eo, on the other hand, the discrete structures due 
to the finite-size effect are all gone, and the susceptibility 
approximates the bulk limit independently of the shape 
and the edge configuration of the flake. Considering eo 
is written as 8000[K]/-\/£[nm] using the graphene's band 
velocity, we find that the room temperature belongs to 
the low temperature regime for a flake of a few nanome- 
ter, while it is in the high temperature regime for a flake 
more than 50nm. 

In the low-temperature regime, the diamagnetic cur- 
rent circulates entirely on the graphene flakes, reflecting 
the absence of characteristic length scale. As the temper- 
ature increases, the current gradually becomes to circu- 
late only near the edge, with the characteristic depth of 
Acdgc = hv /2-nk B T . The local current distribution along 
the cross section perpendicular to the boundary is insen- 
sitive to the global shape of the flake, but significantly 
different between armchair and zigzag edges. 

We predict that the diamagnetism of graphene can be 
possibly observed using the magnetic-field alignment of 
graphene flakes. We estimated the angle distribution at 
various magnetic fields, and found that a strong align- 
ment can be realized in the feasible magnetic field range 
for flakes of S > (10nm) 2 . 
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